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ABSTRACT 



Context. A mysterious X-ray nebula, showing a remarkably linear geometry, was recently discovered close to the Guitar Nebula, the 
bow-shock nebula associated with B2224+65, which is the fastest pulsar known. The nature of this X-ray feature is unknown, and 
even its association with pulsar B2224+65 is unclear. 

Aims. We attempt to develop a self-consistent scenario to explain the complex phenomenology of this object. 

Methods. We assume that the highest energy electrons accelerated at the termination shock escape from the bow shock and diffuse 
into the ambient medium, where they emit synchrotron X-rays. The linear geometry should reflect the plane-parallel geometry of its 
ambient field. 

Results. We estimate the Lorentz factor of the X-ray emitting electrons and the strength of the magnetic field. The former 10 8 ) 
is close to its maximum possible value, while the latter, at = 45 fiG, is higher than typical interstellar values and must have been 
amplified in some way. The magnetic field must also be turbulent to some degree to trap the electrons sufficiently for synchrotron 
X-ray emission to occur effectively. We propose a self-consistent scenario in which, by some streaming instability, the electrons 
themselves generate a turbulent field in which they then diffuse. Some numerical coincidences are explained, and tests are proposed 
to verify our scenario. 

Conclusions. Electron leaking may be common in the majority of pulsar bow-shock nebulae, even though the X-ray nebulosity in 
general is too diffuse to be detectable. 
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1. Introduction 

Pulsar B2224+65 and its bow-shock nebula (the so- 
called"Guitar Nebula") are quite peculiar and extreme objects. 
The pulsar is characterized by an exceptionally high velocity. 
Its proper motion fi = 182masyr~ 1 translates into a transverse 
velocity V* it = l,609(<i/1.86kpc)kms~ 1 , where 1.86kpc is 
the distance estimated from the pulsar dispersion measure 
(according to Cordes & Lazio [2002 ), which implies that it is 
the fastest pulsar known. Otherwise, PSR B2224+65 behaves 
like a standard radio pulsar with a period P = 0.68 s and a 
spin-down power E = 1.20 x 10 33 ergs~', corresponding to a 
characteristic time t = 1.12xl0 6 yr and a pulsar equatorial 
field B* = 2.60 x 10 12 G. The pulsar timing does not present any 
anomaly. Visible in Balmer lines, the nebula displays a peculiar 
shape (as indicated by its nickname), which does not resemble 
in any way a "well-behaved" bow shock. It presents instead a 
conical head, an elongated neck, and a body consisting of a 
main bubble plus probably a smaller one, for a total angular 
length of about 80". Another record of this object is its small 
bow-shock size (0.06" + 0.02"), marginally resolved with the 
Hubble Space Telescope (Chatterjee & Cordes 120021 1. which 
corresponds to a linear size R\, ovl = 1.7 x 10 15 (t//1.86kpc)cm; 
this unusually small size is a natural consequence of the high 
pulsar speed. 

Previous searches for an X-ray counterpart of the nebula pro- 
vided negative or, most confidently, marginal results (Romani 
et al. 119971 ). Only an analysis of Chandra X-ray data by Hui 
& Becker (2007) (also Cordes et al., unpublished manuscript) 



clearly revealed an extended linear feature whose origin apeared 
to coincide with the pulsar position and extended for over 2' 
(which at a putative distance of 1.86 kpc implies a linear size 
Lf ea of over 1 pc), while remaining collimated (with a transverse 
size of about 20") and not having any apparent sign of bending. 
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Kargaltsev & Pavlov (120081) suggested that this feature was 
physically associated with another X-ray source of unknown na- 
ture. However, even though a coincidence cannot be excluded, 
a physical connection between B2224+65 and the X-ray feature 
appears natural. It is supported by the fact that the linear fea- 
ture points almost exactly towards the pulsar and corresponds 
to a hard non-thermal spectrum, in addition to some numerical 
coincidences outlined below. 

Even if association with the pulsar appears likely, the nature 
of this X-ray feature remains unclear. The most puzzling fact 
is that its orientation deviates from the direction of the pulsar 
motion by = 118°, and it is therefore located outside the bow- 
shock region. A scenario involving a ballistic jet does not ap- 
pear viable to be an explanation of this feature, which, although 
protruding for over 1 pc into the ambient medium, remains a re- 
markably linear structure. To excavate dynamically such a long 
path through the ambient medium without experiencing any ap- 
preciable bending, a jet should be highly energetic and extremely 
well collimated. The requirement that the ambient medium ram 
pressure has not significantly affected the direction of this bal- 
listic motion on a scale more than 2,000 times larger than the 
bow-shock size (which is also affected by the same ambient ram 
pressure), places an extremely tight limit on the collimation an- 
gle of one such jet, which we estimate below. 
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We assume that a constant fraction £ of the pulsar spin-down 
power enters an isotropic wind, while a fraction // + fi < 1) 
is channelled into a cylindrical jet of circular cross-section and 
transverse size Wf ea - While the material in the jet travels a dis- 
tance Lf ea , we also assume that the ambient ram pressure deposits 
an extra transverse momentum, equal to a factor of e times the 
original jet momentum (e must be very small, otherwise the jet 
bending would be appreciable). A comparison between the orig- 
inal momentum of the jet and the additional transverse momen- 
tum implies that: 



e/jE/c = WfeaifeaSin(118°)p am bV« 



(1) 



The size of the bow shock is also determined by the balance be- 
tween stellar wind and ambient medium ram pressures, namely 



(2) 



In the above two equations, yu and £ are the relative conversion 
efficiencies. Combining the above equations, we derive: 
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where we have used Lf ea > 2,000 Rbow Even though this re- 
sult can be partially attenuated by assuming an ambient density 
gradient, a low wind efficiency, and/or by invoking some special 
transient in the pulsar energy release (even though no sign of this 
was evident in the pulsar timing), it is theoretically difficult to 
account for such a small Wf ed /Lf e - d ratio, especially considering 
that the measured value of this ratio is about 0.1-0.2. Therefore, 
even though all the above estimates are approximate and some 
assumptions could be refined, the observed value is many orders 
of magnitude above the theoretical upper limit, and therefore the 
hypothesis of a ballistic jet is difficult to pursue. 

In this paper, we propose an alternative scenario to explain 
the linear X-ray feature and its phenomenology, which is based 
on the idea that high-energy electrons may diffuse away from 
the bow-shock region and interact with the ambient medium. 
The basic ideas behind this scenario and its assumptions are pre- 
sented in the next section. Section 3 investigates the evolution of 
the highest energy electrons in the bow-shock region, while the 
interaction of these electrons with the ambient field is consid- 
ered in Sect. 4. Section 5 presents our conclusions and discusses 
the model predictions that should be tested in a near future. 

2. Basic ideas and assumptions 

The principal feature of the scenario proposed here is that the 
highest energy electrons accelerated at the pulsar wind termi- 
nation shock may escape from the bow-shock region and dif- 
fuse through the ambient medium, where they emit synchrotron 
X-rays. These electrons interact directly with the ambient mag- 
netic field, and this interaction will affect both their motion and 
emission. The macroscopic dynamical properties of the ambient 
medium will not be changed; for instance, the orientation of the 
X-ray feature simply reflects the original orientation itself of the 
ambient magnetic field. However, the electrons may play a role 
in the amplification of the ambient magnetic field by creating a 
turbulent component, which may affect their diffusion. The dif- 
fusion coefficient perpendicular to the original orientation of the 
magnetic field remains small, and the cross-field diffusion can be 
neglected. In the following, we verify the internal consistency of 
these assumptions, as well as their consistency with the observed 
phenomenology. 



The pulsar moves with respect to the ambient medium and, 
as a consequence, electrons are always injected in different flux 
tubes. Even in the absence of cross-field diffusion, the nebu- 
lar source has a thickness that depends on the synchrotron life- 
time of the X-ray emitting electrons. The map shown by Hui 
and Becker (120071 1 represents a transverse size of at most 20"; 
while an average transverse profile (taken from Cordes et al., in 
preparation; Romani, private communication) is consistent with 
a sharp (< 2") leading edge and a backward tail of a total thick- 
ness = 18". A fit to this profile by an exponentially decreasing 
law provides an e-fold scalelength = 19". Using this value, the 
synchrotron timescale for the X-ray emitting electrons can be 
estimated to be: 



hit - 



19" 



s 120 yr. 



sin(118°)0.182"yr- 1 
Their lifetime, derived from synchrotron theory, is: 
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to = 24.5 jrVV* 120 ( 

Therefore, the synchrotron timescale as inferred from observa- 
tions corresponds to a magnetic field = 45 /jG (here and in the 
following, we use 2 ke V as a reference energy for the observed 
X-ray photons). With this magnetic field, the Lorentz factor of 
the X-ray emitting electrons is: 



45juG 



3. Physical conditions in the bow-shock region 

The Lorentz factor of the electrons emitting X-rays in the neb- 
ular feature is so large that their gyration radius is comparable 
with the bow-shock size. In fact, if the magnetic field in the head 
of the bow shock is of the order of the equipartition field 



B bow = ^E/cR 2 ^ s> 170£ 1/2 (d/lMkpcT 1 fiG , (7) 

the gyration radius of electrons with Lorentz factor yx is 
# g yr,bow = {m e c 2 1 eB hovl )y x , namely: 
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The fact that /?g yii bow/^bow is close to unity supports the origi- 
nal assumption that high-energy electrons may escape from the 
bow-shock region. The spectrum of the X-ray feature is quite 
hard (as presented by Hui & Becker 2007). This implies that the 
electron energy distribution is dominated by its highest-energy 
part, and that a leakage of the highest-energy electrons would 
affect substantially the entire energy budget of the system. 

Another interesting result is that the electrons produced at 
the termination shock may reach quite high energies. As a di- 
mensional scaling, we evaluate the Lorentz factor of maximally 
accelerated electrons. According to Goldreich & Julian (1969), 
the maximum potential drop (between the pole and the last open 
field line) in an aligned pulsar is: 



Ad) = (aO/c) z flB, = y3£/2c, (9) 
which corresponds to an acceleration of up to a Lorentz factor: 



Tmpd - 



2m e c 2 2m e c 2 



\— =7.2xl0 7 . 
2c 



(10) 
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The processes that lead to the acceleration of high-energy elec- 
trons in pulsar wind nebulae are complex and poorly understood, 
and involve both pulsar wind acceleration and particle accelera- 
tion at the wind termination shock. It appears that 7mpd can be 
used broadly as a reference value: in the Crab Nebula, for in- 
stance, the maximum Lorentz factor of the injected electrons is 
about 10% of 7mpd (De Jager et al. 1 1996k while in this case a 
maximum energy of the injected electrons of the order of 7mpd 
appears to be required. 

4. Physical conditions and processes in the nebular 
feature 

We consider the magnetic field in the nebular feature, whose 
value (45 /iG) was estimated from the transverse profile of the 
X-ray feature. An underlying assumption was that the electrons 
do not diffuse orthogonally in the ambient magnetic field. The 
unavoidable jittering across the flux tube, of the order of the gy- 
ration radius, can be estimated to be: 

* s0 .i3"(-^)(— fM— f ■ (ID 

gy ' \2keV/\45^G/ \1.86kpc/ v 

This value is consistent with the upper limit of 2" to the 
leading edge of the X-ray feature (Cordes et al., unpublished 
manuscript); given the strong dependence on B, a magnetic field 
far lower than our estimate would however imply a R gyl j ea in 
excess of the (measured) leading edge thickness. 

Concerning the origin of this field, one cannot exclude in 
principle the presence of pre-existing inhomogeneities in the am- 
bient medium, as long as they preserve the general plane-parallel 
magnetic-field structure. A filament is clearly visible in the Ho- 
images of the Guitar Nebula, about 100" to the south of and ap- 
proximately parallel to the X-ray feature, and it may be associ- 
ated with a density enhancement (although it is not even known 
whether it is at the same distance of the Guitar Nebula). 

Another more natural possibility is that some turbulent am- 
plification of the originally plane-parallel field has occurred. 
Some degree of turbulence would also help to explain the length 
of the X-ray feature. If electrons with lifetimes of about 120 yr 
(as estimated above) flow along the field lines close to the speed 
of light, they would travel for about 36 pc; the observed length 
of the X-ray feature is, however, only about 1 pc, and it is in- 
sufficient that the electrons become invisible beyond this dis- 
tance. The X-ray luminosity of the nebular feature is extremely 
high, 4.1 x 10 31 (<// 1.86 kpc) 2 erg s" 1 (0.5-10keV band). This 
is 3.5 X 10~ 2 (c//1.86 kpc) 2 times the pulsar spin-down power, 
a high efficiency when compared with those usually measured 
for pulsars (10~ 3 for the pulsar+nebula X-ray emission; Becker 
& Triimper |1997l l, and therefore it is more natural to expect that 
the high-energy electrons lose energy almost completely in the 
X-ray feature. 

To release most of their energy in this distance, the electrons 
cannot flow away freely along the flux tube. A natural alternative 
is that they scatter back and forth along the flux tube: this may be 
possible if, in addition to the plane-parallel large-scale magnetic 
field, there is a turbulent field component. The fact that the linear 
brightness of the nebular feature was observed to decrease "lin- 
early" (Hui & Becker [20071 ) may in fact be reconciled with the 
case of diffusion. In the case of steady injection and diffusion, 
the expected profile resembles a decreasing exponential, and the 
observed profile is consistent with an exponential decrease of 
e-fold scalelength of about 140". In the following, we use this 
value for Lf ea sin t where, for generality, we have also introduced 



the angle t between the direction of the ambient magnetic field 
and the line of sight. Anyway, we do not assume any special field 
orientation in our direction, so that sin t ~ 1 . 

The scattering mean free path along the field lines is: 

/ly * L 2 f JL m a 5" sinr- 2 (d/1.86 kpc) . (12) 



According to the diffusion theory, A|| = ^ gyr with r\ = (<JZ?/B) re 2 , 
where "res" means fluctuations SB of a wavelength that is reso- 
nant with R gyr . Therefore, one can estimate: 




As required for self-consistency, 77 is > 1 and 6B/B < 1. The 
small value for 6B/B refers only to resonant fluctuations and 
may therefore still be consistent with a global magnetic amplifi- 
cation of a few times. 

The derived value for 6B/B is far higher than typical val- 
ues for random fluctuations of the Galactic magnetic field with 
wavelengths of about 10 _3 pc (i.e. the gyration radius of the X- 
ray emitting electrons), which is 6B/B ~ 10 5 — 10 4 . We used 
estimates by De Marco et al. (2007) in the parsec range, and as- 
sumed a Kolmogorov spectrum of fluctuations. A possibility is 
that the flow of relativistic electrons itself generates, by some 
streaming instability, the turbulent component of the magnetic 
field responsible for the field amplification as well as the con- 
finement of the electrons themselves. Without a detailed analy- 
sis of the instabilities, which is beyond the scope of the present 
work, one can verify whether this scenario is energetically con- 
sistent. The power that feeds magnetic fields can be estimated to 
be: 

(B 2 /87r)/? gyr , fea L fea sin r" 1 sin(118 )K,t 

■'•"^"'Issffsv)'"^'"' ,15) 

which is close in value to the X-ray luminosity of the feature, 
4.1 x 10 31 (rf/1.86 kpc) 2 ergs -1 , as if an equipartition between 
electrons and magnetic field has almost been established, on 
timescales shorter than the synchrotron lifetime. 



5. Conclusions 

We have shown that a self-consistent scenario may exist that ex- 
plains the basic attributes of the phenomenology of the X-ray 
feature detected close to the Guitar Nebula. 

Some "numerical coincidences" between measured quanti- 
ties can be readily explained by this scenario. Its assumptions 
could therefore help us to understand other peculiar features in 
the Guitar Nebula, and other less extreme pulsar bow-shock neb- 
ulae, for example the possibility of the highest energy electrons 
escaping from the bow-shock region. 

If the magnetic field within the bow shock is close to equipar- 
tition, the fact that the bow-shock size is similar to the gyra- 
tion radius of electrons with Lorentz factors close to Ympd is 
not just a "numerical coincidence" valid only in the case of the 
Guitar Nebula. The similarity in the two scalelengths is found for 
any pulsar bow-shock nebula, independently of the pulsar spin 
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down-power and velocity. By substituting Eqs. 7 and 10 into the References 
definition of 7?™-, we find that 



m e c 

R gyr = — Tmpd = 

" D bow 



2B h 



1 2c ~ 4 bow ' 



(16) 



The level of proximity of both the bow-shock magnetic field to 
its equipartition value, and the maximum electron Lorentz factor 
to the value ^mpd could produce the observed differences in the 
pulsar bow-shock nebulae. The Guitar Nebula is, however, an 
ideal object in which to observe these effects, because in other 
bow-shock nebulae a similar X-ray feature (if scaled to the bow- 
shock stand-off distance) would be too diffuse to be detectable. 

From the theoretical side, a more detailed investigation is 
required to understand which instabilities could allow an effi- 
cient, turbulent amplification of the magnetic field, as required 
to explain the X-ray feature. It would also be important to deter- 
mine the highest energies of the electrons injected from the wind 
termination shock, a long-standing problem that has never been 
resolved. 

From the observational side, there are some predictions of 
the present model that need to be verified. First of all, the linear 
feature must travel at the same speed as the pulsar: this predic- 
tion could be tested in the next few years because the pulsar, with 
its high proper motion, covers about 1 " in only 5 years. If no mo- 
tion was be detected, a completely different scenario would have 
to be envisaged. Another prediction is the spectral change across 
the feature, with harder spectra on its front side and softer spec- 
tra on its back side (due to the synchrotron-driven evolution of 
the electrons, injected at different times). 

Inverse Compton scattering of electrons with Lorentz factors 
~ 10 s can upscatter CMB photons to the TeV range. Therefore, 
in principle, this object could also be a TeV source, but it would 
be far too faint to be detected by present day Cherenkov tele- 
scopes. At B — 45 /vG, the magnetic field energy density is 
about a factor of 200 higher than the energy density of the 
CMB radiation. Therefore, the total inverse Compton luminosity 
of this feature should be 200 times lower than the synchrotron 
luminosity, namely ^ 2 x 10 29 ergs~' corresponding to a flux 
- 3 x 10~ 16 cm 4 s _l , about three orders of magnitude too faint 
to be detected by MAGIC (5<x detection in 50 hr). Of course, de- 
tecting a TeV source at the position of the feature (its size would 
allow it to be marginally resolved from the pulsar itself) would 
require a revision of the above model. 

Last but not least, the present model does not explain why 
we do not see another X-ray feature on the opposite side of the 
pulsar. In principle, electrons could flow equally well on both 
sides of a magnetic flux tube. Also, no Doppler boosting effect 
can be invoked because, due to scattering, the bulk motion of the 
electrons in the feature is non-relativistic. 

We propose that this asymmetry reflects an asymmetry in the 
pulsar wind itself. Even for an axisymmetric pulsar wind, it is 
sufficient that the symmetry axis of the pulsar wind is not paral- 
lel to the pulsar velocity for such an asymmetry to be produced. 
A detailed numerical modelling is required to compute quantita- 
tively, for different geometries, the resulting level of asymmetry 
in the X-ray source. 
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